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ABSTRACT 

We explore the gas-to-dust mass ratio (Mgas/Md) and the CO luminosity- 
to-Mgas conversion factor (aco) of two well studied galaxies in the GOODS-N 
field, that are expected to have different star forming modes, the starburst GN20 
aX z = 4.05 and the normal star-forming galaxy BzK-21000 at 2 = 1.52. De- 
tailed sampling is available for their Ray leigh- Jeans emission via ground based 
mm interferometry (1.1 — 6.6 mm) along with Herschel PACS and SPIRE data 
that probe the peak of their infrared emission. Using the physically motivated 
Draine & Li (2007) models, as well as a modified black body function, we mea- 
sure the dust mass (Mdust) of the sources and find (2.0]i;Q;g x 10^) M0 for GN20 
and (8.6lo:9 x 10^) Mq for BzK-21000. The addition of mm data reduces the 
uncertainties of the derived Mdust by a factor of ~ 2, allowing the use of the local 
Mgas/Md vs metallicity relation to place constraints on the aco values of the 
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two sources. For GN20 we derive a conversion factor of aco < 1-0 Mq pc~^(K 
km s^^)^^, consistent with that of local ULIRGs, while for BzK-21000 we find a 
considerably higher value, aco ~ 4.0 MqPc~^(K km s~^)~^, in agreement with 
an independent kinematic derivation reported previously. The implied star for- 
mation efficiency is ~ 25 Lq/M© for BzK-21000, a factor of ~ 5 — 10 lower than 
that of GN20. The findings for these two sources support the existence of differ- 
ent disk-like and starburst star-formation modes in distant galaxies, although a 
larger sample is required to draw statistically robust results. 



1. Introduction 

The determination of the conversion factor from CO luminosities to the molecular gas 
mass (Mgas) of a galaxy, (accS = ^gas / -^^co) remains an open issue as there is evidence that 
it varies considerably as a function of metallicity and intensity of the radiation field. Downes 
& Solomon (1998) showed that aco is ^ factor of ~ 6 smaller for local ultra-luminous infrared 
galaxies (ULIRGs, Ljr > lO^^L©) than for local spiral galaxies. A similar picture seems to 
emerge at high redshift, with a fraction of submillimeter galaxies (SMGs) playing the role of 
local starbursts, having lower aco values and exhibiting enhanced star formation efficiencies 
(defined as LiR/Mgas), when compared to normal star- forming galaxies selected by their 
rest-frame UV or optical light (Tacconi et al. 2008, Daddi et al. 2010, Genzel et al. 2010, 
Narayanan et al. 2011). Daddi et al. (2010a,b) used a kinematic analysis of star-forming 
disk galaxies at 2; ~ 1.5 to argue that they have CO conversion factors aco = 3.6 ± 0.8, like 
that in the Milky Way. Instead, Tacconi et al. (2008) and Carilli et al. (2010), placed an 
upper limit of ~ 0.8 on the aco value of high— 2; SMGs. Together with a variety of other 
evidence, these results support the existence of two distinct star formation regimes: a long- 
lasting mode for normal-disks galaxies and a more rapid mode for local starbursts and SMGs 
(Daddi et al. 2008, 2010; Genzel et al. 2010). Nevertheless, excitation biases introduced by 
the use of different molecular lines, along with substantial uncertainties on the aco values, 
raise potential concerns about the role of SMGs in this picture (e.g., Ivison et al. 2011). 

Several studies in the local Universe have tried to tackle this question by measuring the 



^The units of aco, Mopc ^ (K km s ^) ^, are omitted from the text for clarity. 

^Throughout this paper, we use the term "normal" to refer to star-forming galaxies that fall within the 
so-called "main sequence" relation between their star formation rates and stellar masses (e.g., Brinchmann 
et al. 2004; Noeske et al. 2007; Elbaz et al. (2007,2011); Magdis et al. 2010a). Instead, "starbursts" are 
galaxies with substantially elevated SFRs for their stellar masses. 



-3- 



total dust mass of a galaxy (Mdust) and assuming that it is proportional to Mg^s (e.g., Leroy 
et al. 2011). However, determining M^ust is a complex task. In order to break degeneracies 
inherent in current models, a proper characterization of the peak of the spectral energy 
distribution (SED) of the source as well as of the Ray leigh- Jeans emission tail is required. 
With the wide wavelength coverage (70 to 500 fim ) of the Herschel Space Observatory 
{Herschel, Pilbratt et al. 2010), we can now directly observe the peak of the IR emission 
of high— 2; galaxies. Together with ground based mm observations that probe the Rayleigh- 
Jeans emission, this allows us to properly quantify M^ust, and explore possible differences 
between the shape of the SED, the Mgas/Md ratios and the aco values of normal and starburst 
galaxies at high redshift. 

As a test case, we combine HerschelPACS and SPIRE data, from the GOODS- ifersc/ie/ pro- 
gram with (sub)mm observations for two of the best-studied high— 2; galaxies in the Great 
Observatories Origins Deep Survey North field (GOODS-N), the SMG GN20 at z = 4.05, 
and the normal star-forming galaxy, BzK-21000, ai z = 1.521. The uniqueness of these 
sources relies on the detailed sampling of their Rayleigh- Jeans emission via ground based 
mm interferometry (1.1 — 6.6mm). Our aim is to derive Mgas/Mj ratio estimates, investi- 
gate the slope of their SED in the Rayleigh- Jeans tail and put constraints on the aco values. 
Throughout this paper we assume fl^ = 0.3, Hq = 71 km sec~^ Mpc~^, f^A = 0.7 and 
Chabrier IMF. 



2. Sample and observations 

We use deep 100- and 160 /xm PACS and 250-, 350-, 500 ^um SPIRE observations of 
GOODS-N from the GOODS-H program. Details about the observations are given in Elbaz 
et al. (2011). Herschel fluxes are derived from PSF fitting using galfit (Peng et al. 2002). A 
very extensive set of priors is used for 100-, 160- and 250 /xm, including all galaxies detected 
in the ultra-deep Spitzer MIPS 24 /xm imaging, which effectively allow to obtain robust flux 
estimates for relatively isolated sources, even beyond formal confusion limits at 250 /im. For 
350- and 500 /im this approach does not allow accurate measurements due to the increasingly 
large PSFs. Hence wc use a reduced set of priors based on VLA radio detections, resulting 
in flux uncertainties consistent with the confusion noise at these wavelengths. We note that 
GN20, the radio priors include also the nearby GN20.2 objects (both "a" and "b" components; 
Daddi et al. 2009). A detailed description of the flux measurements and Monte Carlo 
derivations of the uncertainties will be presented elsewhere (Daddi et al. in preparation). 

Originally detected at 850 /xm by Pope et al. (2006), GN20 is one of the best studied 
SMGs to date, the most luminous and also one of the most distant {z = 4.055; Daddi et 
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al. 2009a) in the GOODS-N field. It is detected in all Herschel bands apart from 100 jim.. 
Carilli et al. (2010), reported the detection of the C0[1— 0] and C0[2— 1] lines with the Very 
Large Array (VLA), and CO [6— 5] and CO [5— 4] lines with the Plateau de Bure Interfer- 
ometer (PdBI) and the Combined Array for Research in Millimeter Astronomy (CARMA) 
respectively. The source is detected in the Aztec 1.1mm map and continuum emission is 
also measured at 2.2-, 3.3-, and 6.6 mm (Carilli et al. 2011) and at 1.4 GHz with the VLA 
(Morrison et al. 2009). GN20 is identified in all IRAC bands and at 24 //m while it appears 
as a S-band dropout in the ACS-HST image. The stellar mass of the source is 2.3 x 10^^ 
M0 (Daddi et al. 2009a). A compilation of the photometric data is given in Table 1. 

BzK-21000 is a near-IR selected star- forming galaxy, with a spectroscopic redshift, z = 
1.521 (Daddi et al. 2008, 2010a). The source has secure detections in both PACS and in 
the first two SPIRE bands, while at 500 /xm it is only marginally detected. In addition to 
IRAC and MIPS 24 //m data, the source is seen in the 16 //m IRS peak-up image (Tephtz et 
al. 2011). With follow-up VLA and PdBI observations, Dannerbauer et al. (2009), Aravcna 
et al. (2010) and Daddi et al. (2008), have reported the detection of CO[3-2], CO[1-0] 
and C0[2— 1] emission lines respectively. Continuum detections and upper limits are also 
obtained at 1.1-, 2.2- and 3.3 mm (Daddi et al. 2010, Dannerbauer et al. 2009). The stellar 
mass is 7.8 x 10^° Mq (Daddi et al. 2010a). The UV rest-frame morphology of this galaxy, 
the double-peaked CO profile, the large spatial extent of the CO reservoir, and the low 
gas excitation all provide strong evidence that this galaxy is a large, clumpy, rotating disk 
(Daddi et al. 2010a). 



3. Estimating total dust masses 

We employ two methods to derive the dust mass of the galaxies: the physically motivated 
dust models of Draine & Li (2007) (DL07 hereafter), and a more simplistic, but widely used, 
modified black body model (MBB). 

The DL07 models describe the interstellar dust as a mixture of carbonaceous grains 
and amorphous silicate grains. The properties of these grains are parametrized by the PAH 
index, Qpah, defined as the fraction of the dust mass in the form of PAH grains. The majority 
of the dust is supposed to be located in the diffuse ISM, heated by a radiation field with 
a constant intensity C/min- A smaller fraction 7 of the dust is exposed to starlight with 
intensities ranging from Umm to Umax, representing the dust enclosed in photo-dissociation 
regions (PDRs). Following the prescription of DL07, we fit the rest-frame mid-IR to mm 
data points and search for the best-fit model by minimizing the reduced x^- The total dust 
mass is derived from the best-fit model and its uncertainty is estimated by the distribution of 
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Md values that correspond to models with < Xmm + 1 (Avni 1976). Although a different 
grain size distribution would result in different M^nst, we choose to adopt the one prescribed 
by the DL07 models as they can successfully reproduce the IR and submillimeter emission 
for a sample of SINGS galaxies (including both normal and starburst galaxies). 

We also fit the SEDs of our sources with the standard form of a single temperature (Td) 
modified black body, leaving the effective emissivity {/3es) as a free parameter along with 
the dust temperature (Td). For the fit we only consider data points with A^est > 40 /im, to 
avoid emission from very small grains and from the best fit model we can then estimate the 
Afdust from the relation: 

= n ^ ^ R a tT 

where Si, is the observed flux density, is the luminosity distance and Krest = '^oi^o/ Kest)'^ 
is the rest frame dust mass absorption coefficient at the observed wavelength (Li & Draine 
2001). The uncertainty in Mdust is obtained as for the DL07 models. The photometric data, 
along with the best fitting DL07 and MBB models are shown in Fig. 1, while in Table 2 we 
summarize the derived parameters. 

The two methods return Mdust estimates that are in broad agreement. In particular, 
for GN20 we derive Mdust = '^.Otol x 10^ M0(DLO7) and l.bto.t x M0(MBB) while for 
BzK-21000 the corresponding values are Mdust = S.G^gJ x 10* and 7.6t\i x 10*^ M©. Best 
fit MBB models also indicate Td = 33K and 34K for GN20 and BzK-21000 respectively and 
/3efr = 2.1 ± 0.2 and 1.4 ± 0.2. The dependence of Mdust on the other two free parameters 
of the MBB model^, i.e. Td and /3efr, along with the 68% and 99% confidence intervals, are 
shown in Fig. 2. 

To evaluate the significance of adding mm data in the derivation of Mdust, we repeat 
the fitting procedure, this time excluding any data at wavelengths longer than 850 /xm. For 
the DL07 models, the best fit Mdust are unaffected, but the uncertainties increase by a factor 
of ~ 2. Similar results are presented by several studies in the local universe (e.g. Draine et 
al. 2007, Galametz et al. 2011), where they find that in the absence of rest-frame (sub)mm 
data, the derived Mdust estimates and highly uncertain. 



■^These values are derived under the assumption of optical thinness. If we drop this assumption the 
corresponding values are Td =33.8 and 46.3 K for BzK2100 and GN20. 



-6- 



4. Discussion 

We derive a total IR luminosity of Lir = (1.9 ± 0.4) x 10^^ Lq for GN20 and Lir = 
(2.1 ±0.3) X 10^2 for BzK-21000, that correspond to SFRs of ~ 2000 and ~ 210 Mq yr^^ 
and specific star formation rates (sSFR) of ~ 8.6 and ~ 2.6 Gyr~^ Several studies have 
shown that star-forming galaxies at any redshift follow a tight SFR— relation, with outliers 
being starburst galaxies (e.g., Brinchmann et al. 2004, Elbaz et al. 2007, Daddi et al. 2007, 
2009, Magdis et al. 2010a). Our results confirm the sSFR of BzK-21000 is similar to that of 
main sequence galaxies defined in the SFR— M* space at this redshift, while GN20 is located 
in the starburst regime. 

As illustrated in Fig. 1, it appears that GN20 has a steeper slope in the Ray leigh- Jeans 
tail when compared to that of BzK-21000. Despite the large uncertainties, the MBB analysis 
indicates that the PeS values of two galaxies are different at a ~ 2(7 significance level. This 
difference in P^s explains the similar effective derived for the two sources despite the fact 
that the SED of GN20 peaks at - 90 fim while that of BzK-21000 peaks at ~ 100 /xm. The 
best-fitting DL07 models suggest that BzK-21000 has properties similar to those found by 
Draine et al. (2007) for local spirals, with a larger fraction of dust in PAHs (gpAH = 3.9%) 
and a less intense radiation field {Umin = 8) compared to that in GN20 (^pah = 1-12%, 
Umin — 25). Similar results are obtained by Elbaz et al. (2011). 

Several studies have shown that there is a correlation between Mgas/Md and the enrich- 
ment of the ISM of a galaxy. In Fig. 3 (left) we plot Mgas/Md vs metaUicity for a local 
sample studied by Leroy et al. (2011). For consistency we have computed all metaUicities 
on Pettini & Pagel 2004 (PP04) scale. This tight correlation (dispersion of ~ 0.15 dex) can 
be used as a tool to constrain aco- In particular, if we know the metaUicity of a galaxy and 
have measured its Mdust, then we can estimate its Mgas and subsequently the aco value of 
the source, if Lqq is known. In what follows we will attempt to apply this approach for the 
two galaxies in our sample, under the assumption that the observed Mggg/Md— metaUicity 
relation for local galaxies holds at high redshift. 

Having derived estimates for the Mdust values of our sources, we need information on 
their metaUicities, for which we have to rely on indirect indicators. One of these is the M* — 
metaUicity relation of Erb et al. (2006), based on which BzK-21000 has a slightly sub-solar 
metaUicity, Z = 8.65. A similar derivation would follow using the fundamental metaUicity 
relation (FMR) of Mannucci et al. (2010) that relates the SFR and the stellar mass to 
metaUicity. 

For GN20, the situation is somewhat more complicated, as the Erb et al. (2006) relation 
is not sufficiently sampled at the high mass end (> 10^^ M©). Using the FMR, and accounting 
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for the evolution observed for galaxies at 2; > 2.5 (Mannucci et al. 2010), we derive a 
metallicity of Z = 8.8. Another metallicity estimate can be obtained by assuming that 
the huge SFR of GN20 could stem from the final burst of star formation, triggered by a 
major merger that will eventually transform the galaxy into a massive elliptical. Once star 
formation ceases, the mass and metallicity of the resulting galaxy will not change further, 
and one might therefore apply the mass-metallicity relation for present-day elliptical galaxies. 
In this scenario, the metaUicity of GN20 could range between Z — 8.8 and 9.2, although 
a moderately super-solar metallicity is more probable. For our purposes we will adopt a 
metallicity of 8.65 ± 0.2 for BzK-21000 and a whole range of Z = 8.8 - 9.2 for GN20. 

The fit to the local M^s^^/ M^—Z relation (Fig. 3 left), indicates a value of Mg^^/M^ ~ 104 
for BzK-21000 and ~ 75 (35) for GN20 assuming a metaUicity of Z = 8.8(9.2). Santini et 
al. (2010) report similar Mgas/M^ ~ 50 for SMGs. Based on these ratios we first calculate 
Mgas for each galaxy and find 8.9 x 10^° Mq and 1.5 x 10^^ (7.0 x 10^° M©) for BzK- 
21000 and GN20 respectively. Then based on the relation Mgas = aco x -^co estimate 
the corresponding «co factors and find «co = 4.1^2 'y BzK-21000 and in the case of 
GN20 aco = 0.9toi for Z = 8.8 and 0.4t°J for Z = 9.2. The derived values, along with the 
estimates of Leroy et al. (2011) for a local sample, are shown in Fig. 3 (middle). We also 
overplot a sample of local ULIRGs from Solomon et al. (1997), for which we were able to 
compute their metaUicities on the PP04 scale. The quoted uncertainties account both for 
the dispersion of the Mg^s/Ma—Z relation and the uncertainties in Mgas /Ma and Mjust- 

The derived aco values agree with previous independent estimates. In particular Daddi 
et al. (2010), based on the dynamical masses of a sample of 2; ~ 1.5 — 2.0 BzK galaxies 
(including BzK-21000), argued for an average conversion factor of aco = 3.6 for high— 2; star 
forming disks and reported a value of Mga^ = (8.1 ± 1.4) x 10^° M© for BzK-21000 (see also 
Aravena et al. 2010). Furthermore, CariUi et al. (2010), based on the C0[1— 0] transition 
line, estimated the gas mass of GN20 to be Mgas = 1.3 x 10^-*^ x (q;co/0.8) M0, putting 
a coarse upper limit on the conversion factor of aco ~ 0-8 from dynamical constraints. 
This indicates that our aco < 1-0 estimate for GN20 is reasonable. Inverting this line of 
reasoning, placing the sources on the Mgas/Md — metallicity plane of Fig. 3 (left) based on 
the Mgas estimates from the literature, indicates that BzK-21000 is very close to the relation 
defined by a linear regression fit to the local sample. Similarly, GN20 appears to broadly 
follow the local trend. Another way to show that the Mgas/Md— Z relation holds for high-z 
galaxies, is to plot direct observables, without any assumptions for aco- Indeed, in Fig. 3 
(right), we plot L'co/Miust vs metaUicity for our sample and find again that both GN20 and 
BzK21000 follow the local trend. Finally, recent results from Genzel et al. (2011) seem to 
verify our assumption. 
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Despite the substantial uncertainties, the agreement with independent aco estimates 
is reassuring. In the local Universe there is observational and theoretical evidence that 
aco ill starburst galaxies is significantly smaller than that in the Milky Way disk (Downes 
& Solomon 1998; Scovillc ct al. 1997). The derived conversion factor for GN20 is consistent 
with that of local ULIRGs or even lower, while for BzK-21000 the value is considerably 
higher and close to that of local spirals. We stress that without the addition of the mm data, 
the large uncertainties in M^ust would not allow us to derive any meaningful conclusions. 
Furthermore, although the sources appear to have comparable L^i/Lqq ~ 100 Lq (K Km 
s'^ pc^)~^, their SFEs, are considerably different. In particular, for BzK-21000 we derive a 
SEE ~ 25 while for GN20 SEE ~ 100 - 200 (depending on the assumed metallicity).The two 
sources also have comparable M*/Afdust ~ 100. This is also the case for local ULIRGs and 
normal SDSS galaxies (da Cunha ct al. 2010), supporting the idea that the property that 
distinguishes starbursts from normal star-forming galaxies is their enhanced SFR (at fixed 
Mgas, Mdust and M,). 

We conclude by noting that our results, although limited to two sources, are in line 
with previous claims that the star-formation mode of BzK-21000 and other high— 2; galaxies 
in the main sequence of the SFR— plane is different than that of most SMGs, and more 
similar to that of local-disks, despite their very large infrared luminosities and star formation 
rates. Additionally, we confirm the validity of the widely adopted ULIRG-like aco factor 
for SMGs (e.g., Tacconi et al. 2008). We have also demonstrated that the combination of 
Herschel with ground based millimetre data provides a powerful tool to investigate the dust 
and gas properties of high— 2; galaxies. A larger sample is needed to extend this investigation 
and draw statistically robust conclusions. 

We thank Alvio Renzini for useful discussions. GEM acknowledges support from Ox- 
ford University. ED acknowledges funding support from ERC-StG grant UPGAL 240039 
and ANR-08-JCJC- 0008. DR acknowledges support from NASA through a Spitzer Space 
Telescope grant. This work is based on observations made with the Herschel Space Ob- 
servatory, a European Space Agency Gornerstone Mission with significant participation by 
NASA. 
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Table 1: Summary of Herschel and (sub)mm data. 
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(2009); Carilli et al. (2011) 



Table 2: Summary of derived properties of GN20 and BzK-21000 with DL07 and MBB 
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Fig. 1.— Observed SED of GN20 (top) and BzK-21000 (bottom) overlaid with the best fit 
Draine & Li (2007) (DL07) models (red) and the best fit single temperature modified black 
body (blue). The black dashed line is the DL07 model without the stellar component, that 
is depicted with an green dotted-dashed line. The black arrows indicate the 3 a upper limit 
at 6.6 mm and 3.3 mm for GN20 and BzK-21000 respectively. 
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Fig. 2. — Uncertainties on the parameters derived by MBB models for the case of GN20 (left) 
and BzK-21000 (right). The plot shows the equivalent 68% and 99% confidence intervals 
on /3eff plotted against Td, as derived by MC simulations. The region enclosed by the 68% 
confidence level contour is colour coded based on the Mdust that corresponds to each set of 
/3es and Tj. The best fit value is denoted with a solid black cross (see also footnote 3). 
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Fig. 3. — Left: M^^/M^ vs metallicity for a sample of galaxies in the Local Group by Leroy 
et al. (2011) (black squares) and local ULIRGs by Solomon et al. (1997) (orange stars). The 
solid black line the best linear regression fit to Leroy's sample and the green shadowed area 
depicts the dispersion of the correlation. Blue and red circles indicate the position of BzK- 
21000 and GN20 respectively, based on previously published Mgas values that correspond to 
aco = 3.6 for BzK-21000 (Daddi et al. 2010) and aco = 0.8 for GN20 (Daddi et al. 2010a, 
Carilli et al. 2010). GN20 is placed at Z = 8.8 and Z = 9.2 which are the lower and upper 
limits for its metallicity. All metallicities are calculated on PP04 scale. Middle: Constraints 
on aco based on the local Mgas/M^—Z relation shown in the left panel. Colors and symbols 
are the same as in the left panel. Right: L^g/Mdust vs metallicity for the same sample of 
galaxies. 



